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Abstract. Widely distributed lizard species that occur in contrasting habitats generally show distinctive morphological and 
thermal traits that are the results of local adaptations via sexual and natural selection. Here, we assessed the effect of sex 
and seasonality on morphology, thermal ecology, and perch height of Sceloporus grammicus of four populations in central 
Mexico. In all populations, males had larger body sizes, heads, and limbs than females. Females had larger abdominal vol- 
umes than males in two of the four populations. Females and males of the tree-boulder population were smaller than those 
of the cactus and grassland populations. Females of the tree-boulder population were also lighter in weight, with shorter 
forearms, than those of the cactus and grassland populations. Assessing thermal ecology, we found sex-specific differences 
in body temperature and behavioural thermoregulation in the grassland population, and seasonal variation in these ther- 
mal traits for lizards in the cactus, tree-boulder, and grassland populations. Lizards in habitats with closed canopy cover 
had lower body temperatures than lizards in more open habitats. Perch height did not vary between sexes or seasons in 
three of the four populations. We found that the intensity of fecundity advantage varies between female populations, and 
that thermal characteristics of lizards show specific relationships with perch height in each population, which is a signal of 
local adaptation of behavioural thermoregulation. 
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Introduction 


Populations of widely distributed lizard species use differ- 
ent habitat types and therefore experience diverse envi- 
ronmental regimes and selective pressures (e.g., LATTAN- 
ZIO et al. 2020). Habitats vary spatially and temporally in 
their physical structure, providing different gradients of 
abiotic and biotic conditions (SMITH & BALLINGER 2001, 
WILLIAMS & MCBRAYER 2015). Additionally, selective 
pressures, such as sexual and natural selection, may dif- 
fer also between habitat types, and these may then shape 
population-specific, contrasting patterns of male and fe- 
male traits, such as morphology (e.g., HERREL et al. 2001), 
thermal ecology (KIEFER et al. 2005), and microhabitat use 
(e.g., WILLIAMS & MCBRAYER 2015, LEYTE-MANRIQUE et 
al. 2017). 

Lizards occupying habitats with dense canopy cover 
are expected to have small heads, slim and long bodies, 
long tails, and relatively short limbs (HERREL et al. 2001). 


On the other hand, lizards living in habitats with an open 
canopy can be expected to have flat bodies and longer 
hindlimbs and toes (e.g., GARLAND & Losos 1994; HER- 
REL et al. 2001). Furthermore, selective pressures are likely 
to affect the anatomy of male and female lizards different- 
ly, resulting in sexual size dimorphism (SSD; SHINE 1989, 
SCHARF & MERI 2013). For instance, in females, natural 
selection can act through a fecundity advantage, favour- 
ing a greater body size or abdominal cavity volume, which 
may enhance the capacity for producing larger clutches or 
young (e.g., SCHARF & MEIRI 2013). 

The thermal ecology of lizards is influenced directly or 
indirectly by intrinsic and extrinsic factors (HUEY 1991, 
ANDREWS 1998). Although some males and females lizards 
do not differ in body temperature or behavioural ther- 
moregulation (e.g., ADOLPH 1990, HUEY & PIANKA 2007), 
others exhibit intersexual variation in such thermal traits 
(KIEFER et al. 2005). This latter pattern could arise because 
of differences in sex-specific social behaviour (ANDREWS 
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1971) and the degree of SSD, which influences the capac- 
ity of heat retention and gain (i.e., thermal inertia) of each 
sex (CAROTHERS et al. 1998). For instance, territorial male 
lizards may have access to higher ambient temperatures 
than females and thus achieve higher body temperatures 
(CAROTHERS et al. 1998, BEAL et al. 2014). Considering that 
microhabitat temperatures are generally higher in habitats 
with an open canopy than in those with closed canopies 
(VITT et al. 2007), lizards in the former might maintain a 
higher degree of behavioural thermoregulation than liz- 
ards living in the latter because of the greater availability 
of basking sites. In addition, seasonal variation in body 
temperature and the behavioural thermoregulation of liz- 
ards could arise from their seasonal use of different micro- 
habitat types and biotic interactions (HUEY & PIANKA 1977, 
ADOLPH 1990, Losos 2009). 

Lizard microhabitat use frequently differs across and 
within populations, and this variation has been attributed 
to ecomorphological associations (Losos et al. 1994, LAT- 
TANZIO et al. 2020), perch height selection (RAMfREZ-BAU- 
TISTA & BENABIB 2001, ASBURY & ADOLPH 2007), differ- 
ences in thermal environment (ORTEGA et al. 2019), or ter- 
ritoriality (ANDREWS 1971, GARCÍA et al. 2010). Intraspecif- 
ic variation in ecomorphological associations could influ- 
ence the thermal environment as is experienced by lizards 
through perch height selection, and consequently their 
degree of behavioural thermoregulation (ADOLPH 1990, 
LATTANZIO et al. 2020). For instance, some lizard spe- 
cies that inhabit closed-canopy habitats use lower perches 
and achieve lower body temperatures than lizards in more 
open-canopy habitats (HUEY 1974, TERRA et al. 2018). 

Sceloporus grammicus (WIEGMANN, 1828) is the third 
most widely distributed lizard species of the genus Scelopo- 
rus, ranging from southern Texas in the United States to 
northern Oaxaca in Mexico (SITES 1982). Previous studies 
have detected variation in morphology and SSD between 
populations that inhabit different habitat types (HERNAN- 
DEZ-SALINAS et al. 2010, LOZANO et al. 2020). Despite that, 
it is as yet unknown whether sexual differences in abdomi- 
nal volume exist in this species and if that might be the out- 
come of a stronger fecundity advantage on the females of 
some populations than in others. Thermal ecology evalu- 
ations in this species have been carried out in populations 
that inhabit high-altitude environments (LEMos-EsPINAL 
& BALLINGER 1995, ANDREWS et al. 1997, DÍAZ DE LA VE- 
GA-PEREZ et al. 2019) and urban areas (WOOLRICH-PINA 
et al. 2006). Despite using contrasting structural and sea- 
sonal environments, it has been stated that the genus Scelo- 
porus exhibits a conservative pattern in its thermal biology 
(ANDREWS 1998). Thus, studies comparing multiple popu- 
lations that inhabit contrasting habitats would permit sup- 
port for, or the rejection of, this hypothesis. Regarding mi- 
crohabitat use, a previous study found similarities in these 
ecological traits between populations using different habi- 
tat types (LEYTE-MANRIQUE et al. 2017). Nevertheless, it 
is unknown whether there is a relationship between mor- 
phological and thermal characteristics and the selection of 
specific perch heights. 
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Here, we analysed sexual differences in morphology, 
thermal ecology and perch height use across four popula- 
tions of S. grammicus in central Mexico. The four popula- 
tions occur in habitats that are clearly structurally different 
(see Study area, below). For morphology, we predicted that 
within populations, (1) males would have larger body siz- 
es and heads, and longer limbs than females, while (2) fe- 
males would have larger abdominal volumes than males. 
Between populations, we expected that (3) lizards in habi- 
tats with more open canopy cover would have shorter bod- 
ies, as well as longer heads and limbs than their conspe- 
cifics in closed-canopy habitats. Regarding thermal ecol- 
ogy, we expected that (4) males would have higher body 
temperatures and greater levels of behavioural thermoreg- 
ulation than females. Within populations, we predicted 
(5) seasonal variation in body temperature and behaviour- 
al thermoregulation as a result of seasonal variation in mi- 
cro-environmental temperatures. Between populations, we 
predicted that (6) lizards in open-canopy habitats would 
reach higher body temperatures and greater levels of be- 
havioural thermoregulation than lizards in closed-canopy 
habitats. With respect to perch height, we expected (7) that 
males would use higher perches than females, and (8) sea- 
sonal variation in this trait might be associated with the 
variation in the thermal and hydric environment. With- 
in each population, we predicted that (9) morphological 
(body size and length of fore and hind limbs) and thermal 
variables (body, air, and substrate temperatures) would be 
positively correlated with perch height; between popula- 
tions, (10) lizards in open habitats would use higher perch- 
es than lizards in closed-canopy habitats. 


Materials and methods 
Study area 


Fieldwork was carried out in three localities in Hidal- 
go state, Mexico: Zacualtipan (20°32’ N, 98°73’ W, 2017 m 
a.s.l.), Tizayuca (19°51 N, 98°57’ W, 2277 m a.s.l.), and Min- 
eral El Chico (20°10 N, 98°43 W, 2362 m a.s.l.); and one 
in Puebla state, Mexico: Coronango (19°08 N, 98°15 W, 
2183 m a.s.l.). In Zacualtipán (hereafter referred to as the 
‘tree population [TP], Fig. S1A), the mean annual tem- 
perature and total annual precipitation were 14.4°C and 
1325.8 mm, respectively (CONAGUA 2010). The climate 
is humid temperate, with a wet season from June through 
October and a dry season from November through May 
(Garcia 2004, CONAGUA 2010). The vegetation type 
here is a cloud forest dominated by Quercus affinis, Pinus 
patula, and Liquidambar macrophyla (RZEDOWSKI 2006). 
In Tizayuca (hereafter ‘cactus population’ [CP], Fig. S1B), 
the mean annual temperature and total annual precipita- 
tion were 15.0°C and 559.1 mm, respectively (CONAGUA 
2010). The climate in this locale is semi-arid, with a wet sea- 
son from May through October and a dry season from No- 
vember through April (García 2004, CONAGUA 2010). 
The vegetation type is xeric scrubland, with less than 50% 
tree cover, and the ground between shrubs is exposed for 
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most of the year (RZEDOWSKI 2006). This location is mainly 
dominated by species of Opuntia, Yucca, and Agave (RZE- 
DOWSKI 2006). In Mineral El Chico (hereafter ‘tree-boulder 
population [TBP], Fig. SiC), the mean annual temperature 
and total annual precipitation were 14.3°C and 1145.6 mm, 
respectively. The climate at this site is subhumid temperate, 
with a wet season from May through October and dry one 
from November through April (Garcia 2004, CONAGUA 
2010). The vegetation type is pine-oak forest, characterized 
by Abies religiosa, Pinus rudis, and Quercus rugosa (RZE- 
DOWSKI 2006). In Coronango (hereafter ‘grassland popula- 
tior [GP], Fig. S1D), the mean annual temperature and to- 
tal annual precipitation were 16.3°C and 845.7 mm, respec- 
tively (CONAGUA 2010). The climate is subhumid temper- 
ate, with a wet season from May through October and a 
dry season from November through April (Garcia 2004, 
CONAGUA 2010). The vegetation type is grassland, mainly 
dominated by Bothriochloa laguroides, Paspalum dilatatum, 
and Setaria parviflora (FLORES-HUITZIL et al. 2020). All 
study sites consisted of several zones that covered areas be- 
tween 1.0 and 1.5 km”. In terms of relative openness of can- 
opy, we consider TP and TBP as closed-canopy habitats and 
CP and GP as open-canopy habitats (Fig. S1). 


Data collection 


We collected data from all populations from March 2018 
to December 2019. Visits to the study areas were made on 
a bimonthly basis. During each visit, we conducted visu- 
al encounter surveys (LOVICH 2012) from 09:00 through 
18:30 h, ensuring that most of the surveys occurred on sun- 
ny days, when the lizards would have higher activity levels. 
Lizards were captured by hand or after stunning them with 
rubber bands (LovicH 2012). This study was conducted ac- 
cording to the ethics and regulations for animal research 
of the Universidad Autónoma del Estado de Hidalgo and 
Guidelines for the Use of Animals (2012). 

For this study, only adult females and males (= 40 mm 
snout-vent length [SVL]; SITES 1982) were considered. The 
captured lizards were sexed by the presence of enlarged 
postcloacal scales in males and their absence in females 
(SITES 1982). All morphometric measurements were taken 
in the field using a digital calliper (+ 0.01 mm; Mitutoyo), 
and body mass was obtained with a spring scale (+ 0.02 g; 
°*Pesola). We followed LEYTE-MANRIQUE et al. (2017) for 
the recording of all these measurements. Because body 
shape differs between groups of lizards, we modified the 
cylinder formula originally established by SHINE (1992) 
and GOODMAN et al. (2009) to adjust it to lizards of the 
genus Sceloporus. Therefore, abdominal volume (AV) was 
obtained through the formula: V = nr’h, where 'r' is half 
of the average of three cross-sectional measures; 1) width 
at the middle of the body (ie., the widest region of the ab- 
domen), 2) width of the pectoral girdle (region where the 
forelimbs are inserted), and 3) width of the pelvic girdle (re- 
gion where the hindlimbs are inserted); and ‘h' is the length 
of the abdomen (measured from the neck to the anterior 


margin of the cloaca). We decided to measure ‘h' from the 
neck region because embryos were found in this region 
in previous dissections of females of Sceloporus grammi- 
cus (A. RAM{REZ-BAUTISTA pers. comm.). These variables 
were measured with the same digital calliper (+ 0.01 mm; 
Mitutoyo) as the other morphometric characters. 

Body temperature (T,) was recorded within 15 seconds 
after the capture ofa lizard, using a Miller- Weber rapid reg- 
istering thermometer (+ 0.2°C; VITT et al. 2008). To avoid 
heat gain from the collector’s hand, lizards were held by the 
head while T, was measured (Virr et al. 2003). We then 
recorded the temperature (T) in the substrate where that 
lizards had been found and the air temperature (T) in the 
shade at 1m above the ground with the same Miller-Weber 
thermometer (RAM{REZ-BAUTISTA & BENABIB 2001, VITT 
et al. 2008). Additionally, we measured perch height with 
a tape measure (+ 0.1 cm) as the vertical distance from the 
ground to the point where the lizard had first been seen 
(RAM{REZ-BAUTISTA & BENABIB 2001). 

To avoid pseudoreplication, we marked each lizard in 
two ways: painting a small temporary mark of nail polish 
on the dorsum at the hindlimb level, and permanent toe- 
clipping according to TINKLE (1967). After measuring mor- 
phological and thermal characteristics, we released each 
lizard at the point where it was first observed. We collected 
data from 336 adult lizards, of which 50 were females and 
23 were males from the TP, 45 females and 28 males from 
the CP, 36 females and 20 males from the TBP, and 69 fe- 
males and 65 males from the GP. From these samples, all 
lizards that had required more than 2 minutes from first 
sighting to capture were excluded from the thermal ecolo- 
gy analyses. Lizards found on the ground were not includ- 
ed in the perch height analyses. 


Data analyses 


To remove the effect of size (SVL), we regressed all mor- 
phological variables against SVL and obtained residuals for 
males and females in each population. Then we used these 
residuals in subsequent analyses as SVL-adjusted variables. 
To assess whether there is SSD in each population, we em- 
ployed the two-tailed Student t-test to compare morpho- 
logical variables between sexes (ZAR 2010). To detect in- 
terpopulational differences in morphological variables 
in females and males, we applied a one-way multivariate 
nonparametric analysis (PERMANOVA) using the Gow- 
er similarity distance, with 9999 random permutations, 
and Bonferroni probability adjustment (ANDERSON 2001, 
QUINN & KEOUGH 2002). 

To quantify the degree of behavioural thermoregula- 
tion (AT), we calculated the absolute values of the differ- 
ence between T, and environmental temperatures (T, and 
T), where AT, = (T, -T,) and AT, = (T, -T) (VRCIBRAD- 
IC & ROCHA 1998, KIEFER et al. 2007, MAIA-CARNEIRO & 
ROCHA 2013, 2020). Higher absolute values of AT, and AT, 
indicate higher degrees of thermoregulation (active ther- 
moregulation) in relation to environmental temperatures 
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(KIEFER et al. 2007, MAIA-CARNEIRO & ROCHA 2020). 
Additionally, the degree of behavioural thermoregula- 
tion was evaluated by calculating the percentage of nega- 
tive values of AT, and AT, which indicates the frequency 
at which body temperatures were lower than the environ- 
mental temperatures (KIEFER et al. 2007, MAIA-CARNEIRO 
& ROCHA 2013, 2020). Low percentages of negative values 
would indicate a lower level of behavioural thermoregula- 
tion, whereas higher percentages suggest a higher degree 
of active thermoregulation (MAIA-CARNEIRO & ROCHA 
2013, 2020). The interpretation of the degree of behavioural 
thermoregulation through AT and percentage of negative 
values relies on knowing which are the main behavioural 
mechanisms for thermoregulation and how they influence 
the difference between T, and environmental temperatures 
(T, and T) (Mata-CARNEIRO & ROCHA 2020). 

To evaluate the effects of seasonality and sex on T,, envi- 
ronmental temperatures (T, and T), and behavioural ther- 
moregulation (AT, and AT) within each population, we 
performed a two-way PERMANOVA on Gower distances 
calculated from raw data with 9999 random permutations 
and a multiple nonlinear regression with season and sex 
as factors (ANDERSON 2001, QUINN & KEOUGH 2002). To 
identify the relationships between T, and environmental 
temperatures, and Ty and SVL, we carried out Spearman 
correlation analyses between these variables (ZAR 2010). 
Because thermal variables did not differ between sexes (in 
three out of four populations) but did differ between sea- 
sons in some populations (see Thermal ecology, below), we 
used a two-way PERMANOVA with similar characteris- 
tics and a multiple nonlinear regression with season and 
population as factors (ANDERSON 2001, QUINN & KEOUGH 
2002). 

To assess the effects of seasonality and sex on perch 
height within each population, we used a two-way analysis 
of variance (ANOVA) with season and sex as independent 
categorical variables, and perch height as the dependent 
variable (ZAR 2010). To identify the relationships between 
morphological and thermal variables and perch height, 
we carried out Pearson correlation analyses between these 
variables with Bonferroni probability adjustment (ZAR 
2010). Because perch height did not differ between sexes 
nor seasons in three out of the four populations (see Perch 
height below), we carried out a one-way ANOVA to es- 
tablish interpopulational differences in perch height (Zar 
2010). Our results are presented as means + SE unless indi- 
cated otherwise. We considered differences with an a value 
of 0.05 to be significant. Statistical analyses were computed 
in Past 4.0 (HAMMER et al. 2001) and Statistica 10.0 (Stat- 
soft, Inc., Tulsa, OK, USA). 


Results 
Morphology 
In the TP, males had larger heads and longer femurs (FL) 


and tibias (TL) than females (Supplementary Table S1). 
However, females had greater abdominal volumes than 
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males (Supplementary Table S1). In the CP and TBP, males 
were larger than females in all traits except body mass and 
abdominal volume (Supplementary Table S1). In the GP, 
males had a larger SVL, larger and wider heads, and longer 
FL and TL than females (Supplementary Table S1); females 
had greater abdominal volumes than males (Supplemen- 
tary Table S1). The PERMANOVA showed up differences 
between female populations (F, 96 = 4-708, P < 0. 001). Fe- 
males of the TBP were smaller and had shorter forearm 
lengths (FOL) than females of the TP (SVL: P < 0.001, 
FOL: P = 0.024), CP (SVL: P < 0.01, FOL: P < 0.05) and GP 
(SVL: P < 0.001, FOL: P < 0.001, Supplementary Table S1). 
Females of the TBP were lighter in weight than those of 
the TP (P < 0.01) and GP (P < 0.001). Females of the CP 
had narrower heads than females of the TP (P < 0.001) and 
GP (P < 0.01; Supplementary Table S1). The PERMANOVA 
also revealed differences between male populations (F, |. = 
2.196, P < 0.01). Males of the TBP were smaller than those 
of the CP (P < 0.01) and GP (P < 0.001). Males of the CP 
had narrower heads than those of the TP (P < 0.05) and GP 
(P < 0.01, Supplementary Table S1). 


Thermal ecology 


In the TP, neither the effect of season (F, 
or sex (F, œ =07LP= 0.18) nor the interaction between 
them (F, „ = —20.75, P = 0.98) on the thermal variables 
(Table 1) ‘turned out to be significant. In both seasons, the 
percentage of negative values of AT, were female-biased 
(Fig. 1A), whereas those of AT, were male-biased (Fig. 2A). 
The T, of both sexes was positively correlated with T, (r, = 
0.57, P < 0.001, n = 73) and T (r, = 0.65, P < 0.001, n = 73) 
but not significantly correlated with SVL (r, = 0.02, P = 
0.89, n = 73). Compared to the TP, the CP experienced a 
significant effect of season (F, ,, = 3.60, P < 0.001), but not 
sex (F_,, = 0.39, P = 0.48), on thermal variables, and the 
interaction between them again was not significant (F, , 
= —14.46, P = 0.85). Mean T, T, T, and AT, varied inter- 
seasonally but not intersexually (see Tables 1 and S2); liz- 
ards had higher values of T, T, and T, in the dry than in 
the wet season, whilst they had higher values of AT, in the 
wet than in the dry season (Table 1). However, AT, did not 
differ seasonally or sexually (Tables 1 and S2). Percentages 
of negative values of AT, were male-biased in both seasons 
(Fig. 1B). On the other hand, those of AT, were female-bi- 
ased and male-biased in the wet season and dry season, 
respectively (Fig. 2B). The T, of both sexes was positively 
correlated with T, (r,= 0.66, P < 0.001, n = 66) and T (r = 
0.73, P < 0.001, n = 66), but not significantly correlated 
with SVL (r, = —0.08, P = 0.53, n = 66). 
In the TBP, a significant effect of season (F, «= 0.89, 
P < 0.05), but not sex (F, 6 = LOL P= 0.06) or the in- 
teraction between them (E = 79-64 P = 0. 93), on the 
thermal variables was found (Table 1). The T, was higher 
in females than in males in the dry season, whereas the 
AT. of both sexes was higher in the wet season than in the 
dry season (Tables 1 and S3). The percentages of negative 
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Table 1. Summary of thermal characteristics of female and male Sceloporus grammicus of four populations in central Mexico sorted 
by season. Ranges for each group are shown below mean + SE. All variables are given in °C. Abbreviations: TP, ‘tree population; 
CP, ‘cactus population’; TBP, ‘tree-boulder population; GP, ‘grassland population’; n, sample size; T, body temperature; T, air tem- 
perature; T, substrate temperature; AT, absolute difference between body and air temperatures; AT, absolute difference between body 
and substrate temperatures. 


T, T, T, AT, AT, n 

TP 

Wet season 

Females 28.58 + 1.43 24.22 + 0.99 24.18 + 0.82 4.88 + 1.08 4.64 + 0.90 10 
21.80-34.40 20.20-28.20 20.20-27.80 0.00-9.40 0.60-9.00 

Males 29.13 + 1.59 24.35 + 1.21 24.08 + 1.05 5.03 + 1.06 5.05 + 0.80 8 
20.20-33.00 16.80-27.40 17.80-26.80 1.00-9.40 1.00-8.00 

Dry Season 

Females 27.02 + 1.02 22.93 + 0.82 22.84 + 0.81 4.45 + 0.49 4.52 + 0.58 40 
14.00-33.60 12.60-32.00 13.60-31.60 0.40-13.00 0.00-14.00 

Males 29.55 + 1.42 25.95 + 1.45 26.21 + 1.30 4.53 + 0.80 3.95 + 0.62 15 
15.40-35.00 12.00-32.00 14.00-34.60 0.60-9.40 1.00-8.20 

CP 

Wet season 

Females 28.27 + 0.92 24.57 + 0.61 25.13 + 0.66 4.53 + 0.60 3.98 + 0.47 26 
18.80-34.80 19.80-32.80 18.60-31.80 0.00-11.20 0.20-9.80 

Males 29.72 + 1.10 24.85 + 0.90 25.31 + 0.90 5.54 + 0.70 5.09 + 0.58 21 
19.20-38.00 17.20-36.00 16.80-34.60 0.00-11.40 0.80-10.80 

Dry Season 

Females 32.90 + 0.50 30.47 + 0.95 30.21 + 0.79 3.26 + 0.56 3.31 + 0.56 14 
29.80-36.40 23.80-36.00 25.00-35.60 0.80-8.60 0.00-7.40 

Males 31.08 + 2.05 30.68 + 1.62 29.16 + 1.69 2.00 + 0.74 2.96 + 0.77 5 
26.40-36.00 26.80-34.00 25.00-33.00 0.20-3.60 1.40-5.60 

TBP 

Wet season 

Females 29.39 + 0.64 24.24 + 0.65 25.24 + 0.85 5.86 + 0.67 5.52 + 0.61 17 
24.00-34.00 20.00-30.00 18.80-31.20 0.00-10.20 1.20-11.40 

Males 30.13 + 0.98 24.50 + 0.84 25.10 0.82 5.63 + 1.12 5.03 + 1.32 6 
26.20-33.20 21.60-27.20 22.80-28.00 2.20-8.80 1.00-8.60 

Dry Season 

Females 30.55 + 0.92 26.39 + 0.98 27.48 + 0.78 4.56 + 0.95 3.45 + 0.57 17 
21.00-34.80 20.00-32.40 21.00-32.20 0.20-13.00 0.80-8.40 

Males 27.43 + 1.22 23.41 + 0.90 23.47 + 1.20 4.01 + 0.98 4.24 + 0.96 14 
21.00-34.80 19.00-29.40 18.40-32.20 1.00-14.40 0.20-14.00 

GP 

Wet season 

Females 32.34 + 0.42 26.05 + 0.51 26.81 + 0.49 6.34 + 0.56 5.68 + 0.52 37 
25.80-35.60 20.40-31.00 22.40-34.00 0.00-13.20 0.40-11.60 

Males 33.48 + 0.35 26.21 + 0.66 27.78 + 0.55 7.29 + 0.75 5.87 + 0.67 34 
28.80-37.00 20.00-34.60 21.20-34.40 0.00-15.00 0.00-14.40 

Dry Season 

Females 30.19 + 0.55 25.34 + 0.59 25.53 + 0.54 5.23 + 0.46 4.76 + 0.43 32 
22.00-35.40 19.00-30.40 18.60-30.40 0.60-10.60 0.40-12.00 

Males 31.45 + 0.61 25.38 + 0.62 25.90 + 0.65 6.48 + 0.58 5.75 + 0.60 31 
23.20-37.40 20.00-31.60 20.20-33.20 1.80-11.80 0.60-12.80 
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Figure 1. Percentages (%) of negative values of AT, for female (white bars) and male (gray bars) Sceloporus grammicus of four popula- 
tions in central Mexico sorted by season. Populations: ‘tree population’ (TP, A), ‘cactus population’ (CP, B), ‘tree-boulder population 


(TBP, C), and ‘grassland population (GL, D). 
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Figure 2. Percentages (%) of negative values of AT, for female (white bars) and male (gray bars) Sceloporus grammicus of four popula- 
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tions in central Mexico sorted by season. Populations: as in Figure 1. 
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values of AT were female-biased in the wet season, whilst 
they were male-biased in the dry season (Fig. 1C). Con- 
versely, those of AT. were female-biased in both seasons 
(Fig. 2C). The T, of both sexes was positively correlated 
with T, (r, = ne. P < 0.01, n = 54) and T, (r = 0.48, P < 
0.001, n = 54), but not significantly correlated with SVL 
(r, = —0.01, P = 0.95, n = 54). In the GP, the effect of sea- 
son was significant on thermal variables (F, yo = 5-08, P < 
0.01), but neither on sex (F, go = 1-82, P = 0.15) nor the in- 
teraction between them (F, wo = 72-06, P = 0. 62). Howev- 
er, T, and T, varied between seasons, but only T, differed 
between sexes (Supplementary Table $4). Mean E and T, 
were higher in the wet season than in the dry season, and 
in both seasons males had a higher T, than females (Ta- 
ble 1). In both seasons, the percentages of negative values 
of AT, were slightly male-biased (Fig. 1D). Those of AT, 
were, however, male-biased in the wet season, but female- 
biased in the dry season (Fig. 2D). As in the other popula- 
tions, the T, of both sexes was positively correlated with T, 
(r = 0.23, P < 0.01, n = 134) and T, (r,= 0.31, P < 0.001, n= 
134), but not significantly correlated with SVL (r, = 0.08, 
P = 0.36, n = 134). 

Highly significant effects of population (F, „= 3.70, P < 
0.001) and season (F, ju = 2-96, P <0. 001) on thermal vari- 
ables were detected in all populations, but the interaction 
between them was not significant (F, su = -18.77 P = 0.45). 
All but one thermal variables varied i interpopulationally, 
with only T, indicating a trend towards significance (Sup- 
plementary Table S5). Mean AT, and AT, differed intersea- 
sonally as well (Supplementary Table S5). Lizards from the 


GP had a higher T, than lizards in the other populations 
(Supplementary Table S6, Fig. 3A). Environmental tem- 
peratures (T, and T) were higher in the GP and CP than 
in the TP (Supplementary Table S6, Figs 3B-C). In both 
seasons, mean AT, was higher in the GP than in the CP 
(Supplementary Table S6). Across populations, mean AT, 
was higher in the GP than in the other populations (Sup- 
plementary Table S6, Fig. 3D). Mean AT. was higher in the 
GP compared with CP (Fig. 3E). Considering all popula- 
tions, both thermoregulation indices were higher in the 
wet (AT: 5.91 + 0.28; AT: 5.23 + 0.24; for both variables 
n = 159) than in the dry season (AT; 4.78 + 0.25; AT; 4.46 
+ 0.23; for both variables n = 168; P < 0.01 for AT, and P < 
0.05 for AT). The highest percentages of negative values of 
AT, were found in the TP and the CP during the wet sea- 
son and dry season, respectively (Supplementary Table S6). 
Considering both seasons, CP lizards showed the highest 
percentages of negative values of AT, (Supplementary Ta- 
ble S6). The lizards of TBP and CP had the highest percent- 
ages of negative values of AT, in the wet season and dry sea- 
son, respectively (Supplementary Table S6). When pooling 
both seasons, TBP lizards showed the highest percentages 
of negative values of AT, (Supplementary Table S6). 


Perch height 


Mean perch height did not differ between sexes and sea- 
sons in the TP (season effect: F, 0.41, P = 0.11; sex effect: 
F = 2.71, P = 0.11; interaction: F, go Z 0-96, P = 0.33), CP 
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Figure 3. Summary of the variation in body (A), air (B), and substrate (C) temperatures, differences between body temperature and air 
temperature (AT , D) and differences between body temperature and substrate temperature (AT , E) between populations of Sceloporus 
grammicus. Populations: as in Figure 1. Asterisks between bars indicate significant differences between groups (*P < 0.05, **P < 0.01, 


*** P < 0.001). Error bars represent a 95% confidence intervals. 
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Table 2. Perch height (cm) of female and male Sceloporus gram- 
micus of four populations in central Mexico sorted by season in 
central Mexico. Ranges for each group are shown below mean + 
SE. Abbreviations are as in Table 1. 


Females n Males n 

TP 

Wet season 170.00 + 16.81 6 104.16 + 26.30 7 
115.00-210.00 37.00-209.00 

Dry season 161:594:12:88 35 144.83 + 27.55 12 
16.00-320.00 19.00-290.00 

CP 

Wet season 87.85 + 15.27 17 87.56 + 14.88 15 
25.60-204.50 15.50-198.10 

Dry season 115.56 + 15.42 10 109.96 + 51.68 5 
43.10-190.00 30.00-300.00 

TBP 

Wet season 124.79 + 19.33 13 154.67 + 13.37 6 
23.60-220.00 110.00-192.00 

Dry season 164.79 + 14.60 15 121.41 + 17.86 11 
77.00-256.00 40.50-224.00 

GP 

Wet season 90.95 + 8.06 36 87.15 + 8.58 32 
28.50-209.10 10.00-223.00 

Dry season 128.88 + 11.68 30 93.55 + 10.55 24 
20.00-323.00 32.00-199.50 

(season effect: F _ = 1.37, P = 0.25; sex effect: F | = 0.02, 


P = 0.89; interaction: F , = 0.02, P = 0.90), and TBP (sea- 


son effect: F. a = 0-03, 23 0.86; sex effect: E, a = 0-13, P= 
0.72; interaction: F, „ = 3-82, P = 0.06) (Table 2). In the TP, 

the snout-vent, forearm, and femur lengths of lizards did 
not correlate with perch height, but tibia length was nega- 
tively correlated with perch height (Table 3). None of the 
thermal variables were correlated with perch height in liz- 
ards of this population (Table 3). Similarly, in the CP and 
TBP, the morphological variables of the lizards did not cor- 
relate with perch height (Table 3). Regarding the relation- 
ship with thermal variables in the CP, all of them were posi- 
tively correlated with perch height (Table 3). On the other 
hand, the T, of lizards in the TBP was positively correlat- 
ed with perch height, whilst T, and T, were not (Table 3). 
In the GP, females used higher perches in the dry than in 
the wet season (E. as = 5-17 P < 0.05), and they used mar- 
ginally higher perches than males (F, „= 4.03, P < 0.05), 

but the interaction between both factors was not significant 
(F, „s = 2-61, P = 0.11) (Table 2). Morphological variables of 
the GP lizards did not correlate (Table 3), but all thermal 
variables were negatively correlated with perch height (Ta- 
ble 3). Additionally, perch height differed between popula- 
tions (F wo = 13-22, P < 0.001). Lizards of the TP (152.38 + 
10.04 cm, n = 60) and TBP (141.28 + 9.01 cm, n = 45) used 
higher perches than lizards in the CP (96.01 + 9.42 cm, n= 
47) and in the GP (99.79 + 5.00 cm, n = 122) (see also Fig. 4). 
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Table 3. Pearson correlation coefficients (r) between perch height 
and morphological and thermal variables in four populations of 
Sceloporus grammicus in central Mexico. Morphological variables 
are SVL, FOL, FL, and TL; all are in mm. Thermal variables are 
Tp To Ts all are in °C. Abbreviations for morphological charac- 
teristics: SVL, snout-vent length; FOL, forearm length; FL, fe- 
mur length; TL, tibia length. Abbreviations for populations and 
thermal characteristics follow Table 1. n, sample size. *P < 0.05, 
*P < 0.01. 


TP CP TBP GP 
Trait / n 60 47 45 122 
SVL 0.08 0.19 0.22 0.08 
FOL -0.14 -0.24 0.04 -0.09 
FL -0.04 -0.01 -0.19 -0.13 
TL -0.26 * -0.20 0.01 -0.11 
T, 0.03 0.40 ** 0.36 * -0.26 ** 
T, -0.07 0.29 * 0.09 -0.19 * 
T, -0.13 0.35 * 0.29 -0.27 ** 
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Figure 4. Variation in perch height between four populations of 
Sceloporus grammicus. Abbreviations follow Figure 1. Asterisks 
between bars indicate significant differences between groups ac- 
cording to Bonferroni post hoc tests (**P < 0.01, ***P < 0.001). 
Error bars represent a 95% confidence intervals. 


Discussion 
Morphology 


Sexual size dimorphism (SSD) was apparent in all lizard 
populations examined and this result fully supports our 
first prediction, but the second prediction only partially. 
Despite differences in habitat use, our analyses indicate 
that the direction of sexual selection in males is similar 
in all populations. However, the intensity of the fecundity 
advantage might be greater in females of the TP and GP 
than those in the CP and TBP. In the Asian lacertid liz- 
ard Takydromus septentrionalis, female populations differ 


Life history variation of Sceloporus grammicus 


in abdominal length, i.e., those that have shorter abdomens 
produce eggs with lower egg masses (Du et al. 2005). In ad- 
dition, body size and reproductive traits in the Asian aga- 
mid Phrynocephalus przewalskii are related to food avail- 
ability in that larger females produce more and larger eggs 
in localities with higher food availability than those where 
less food is available (ZENG et al. 2013). Hence, our results 
might be indirect evidence of a greater food availability in 
the TP and the GP, which potentially promotes female re- 
productive trait variation in S. grammicus; however, this 
assumption still remains to be tested. Despite that, our 
results about intersexual variation in abdominal volume 
should be taken with certain caution, given that, in some 
cases, we could not clearly identify (by abdominal palpa- 
tion) the reproductive status of females with greater ab- 
dominal volumes. 

Males and females of the TBP were smaller, had short- 
er forearm lengths and were lighter in weight than those 
from the other populations. These results partially support 
our third prediction and are consistent with the findings of 
previous studies of S. grammicus (LOZANO et al. 2020) and 
other lizard species (HERREL et al. 2002, HOPKINS & ToL- 
LEY 2011, LATTANZIO et al. 2020), i.e., that lizards in open 
habitats have longer forearms, which enhances their per- 
formance in climbing on irregular and three-dimensional 
substrates (MELVILLE & SWAIN 2000) and moving about in 
sparsely vegetated areas (HOPKINS & TOLLEY 2011). Large- 
bodied and long-limbed females of the TP, CP, and GP 
would be able to have improved sprint speeds, enhancing 
their predation success on specific prey types or escaping 
from potential predators (MELVILLE & SWAIN 2000, Losos 
et al. 2006). Conversely, short-limbed females of the TBP 
probably would not run faster, but could be expected to 
leap farther than their counterparts in the other popula- 
tions. The larger body sizes of males in the open-canopy 
habitats (CP and GP) might be the outcome of the dryness 
and seasonality of these environments, which promotes a 
high intensity of competition for resources (e.g., females, 
shelters, and basking sites) among males, as occurs in 
Anolis lizards (Losos 2009). 


Thermal ecology 


The results of sexual variation in T, of the GP partially 
support our fourth and fifth predictions. Two contrast- 
ing sexual patterns were found in T, of S. grammicus: one 
in which males attain higher T, than females (ANDREWS 
et al. 1997, WOOLRICH-PINa et al. 2006) and another in 
which there is no substantial intersexual variation in T, 
(LEMOS-EsPINAL & BALLINGER 1995, ANDREWS et al. 1997, 
Diaz DE LA VEGA-PEREz et al. 2019). Our results from the 
GP fit the first pattern, whereas those of the TP, CP, and 
TBP fit the second pattern. Although lizards within the 
same population are generally exposed to similar thermal 
environments, they will often select specific microclimat- 
ic conditions. For example, territorial males of S. jarrovi 
have better or more access than females to microhabitats 


with higher thermal quality for thermoregulation, thus 
this circumstance contributes to intersexual differences 
in T, (BEAL et al. 2014). The similarity in T, between the 
sexes in the other three populations concurs with the fact 
that females had higher percentages of negative values 
of AT, and AT, than males, suggesting that they actively 
thermoregulate in a manner similar to males (VRCIBRAD- 
IC & ROCHA 1998, KIEFER et al. 2007, MAIA-CARNEIRO & 
ROCHA 2013). 

In all populations, the T, of lizards was positively cor- 
related with T, and T, but the relationship was always 
stronger with T, Additionally, the correlation coefficients 
of both relationships were higher in the TP and CP than 
in the TBP and GP. Our results confirm that the heat gain 
strategy of S. grammicus consists of both thigmothermy 
and heliothermy, as occurs in populations of urban areas 
(WooLRIcH-PINa et al. 2006) and of high-altitude habitats 
(Diaz DE LA VEGA-PEREZ et al. 2019). During the coolest 
periods of the day, many lizard species assume body pos- 
tures that maximize heat gain (either from thigmothermy 
or heliothermy), whereas in the hottest periods, they adopt 
postures that favour heat loss (GILLIS 1991, MCCONNACHIE 
et al. 2009). Although we did not directly evaluate body 
posture and sun exposure of lizards, our study may support 
the idea that the thigmothermic strategy of S. grammicus 
may vary between populations inhabiting contrasting hab- 
itats, as has been noted in other lizard species (GILLIS 1991, 
McConnachlE et al. 2009). 

The results of interpopulational variation in T, envi- 
ronmental temperatures (T, and T) and behavioural ther- 
moregulation (AT, and AT) partially support our sixth ex- 
pectation and are consistent with previous studies of differ- 
ent lizard genera, which have found that lizards from open- 
canopy habitats are generally exposed to higher environ- 
mental temperatures and thus achieve higher body temper- 
atures than conspecifics in closed-canopy habitats (HUEY 
1974, WILLIAMS & McCBRAYER 2015, TERRA et al. 2018). On 
the other hand, this intraspecific variation contrasts with 
the conservative pattern of thermal ecology found in stud- 
ies of Sceloporus and other lizard genera, probably reflect- 
ing an ancestral condition (GILLIS 1991, ANDREWS 1998). 
This conservative pattern in thermal ecology across popu- 
lations of the same species has been observed even though 
they are found in differently structured habitats and ex- 
perience contrasting thermal environments (ANDREWS 
1998, SALES & FREIRE 2019). In sum, we demonstrated that 
S. grammicus is relatively plastic in its thermal traits, and 
the way it behaviourally thermoregulates will depend on 
the structure of the habitat it uses. 

Throughout our comparative approach, we found that 
most of the morphological and thermal characteristics are 
subject to conspicuous intraspecific variation in S. gram- 
micus. Previous studies in Anolis and Sceloporus lizards 
have uncovered different patterns across species and popu- 
lations: 1) thermal biology exhibits more divergence than 
does morphology, 2) diversification of thermal biology of- 
ten follows diversification in morphology, and 3) thermal 
biology is similar but morphology diverges conspicuously 
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(GUNDERSON et al. 2018, MUNoz & Losos 2018, LATTAN- 
ZIO et al. 2020). In this manner, multiple contrasting se- 
lective pressures (natural and sexual) can operate simulta- 
neously, promoting patterns of both similarity and diver- 
gence, and behaviour may be the bridge that links these 
outcomes (MUNOz & Losos 2018). Despite being larger in 
many morphological characteristics, males did not nec- 
essarily attain a higher T, than females, except in the GP. 
In fact, percentages of negative values of AT, and AT, indi- 
cated that the behavioural thermoregulation of females is 
sometimes similar to, or even as precise as, the behaviour- 
al thermoregulation of males (SALES & FREIRE 2019). This 
is particularly relevant considering that the T, of females 
will enhance offspring fitness by accelerating embryogen- 
esis and thus parturition (SHINE & DOWNES 1999). Thereby, 
the ability of S. grammicus to take advantage of its available 
thermal environments depends on the way it behaviourally 
makes the most of its morphology, which in turn affects the 
preferred type and height of its perch. 


Perch height 


The results of the effects of sex and season on perch height 
did not confirm our seventh prediction, but weakly sup- 
ported our eighth prediction. Both patterns are consist- 
ent with those found in S. melanorhinus, in which females 
use higher perches than males and there is no seasonal 
variation in perch height (Garcia et al. 2010). Hence, the 
females and males of S. grammicus in our study used the 
same perch height throughout the year for foraging, so- 
cial interactions, and in some instances for thermoregula- 
tion (see next paragraph). The seasonal variation in perch 
height by females of the GP could be linked to the sea- 
sonal foraging rate, avoidance of certain social interac- 
tions or potential predators, or even to remain close to 
specific refuges, as occurs in other lizard species (GARCÍA 
et al. 2010). Considering the unique relationship between 
perch height and T, found in each studied population, it 
is possible that each of these aspects is enhanced in fe- 
males and males if they attain or maintain a higher and/ 
or stable T 

In all populations, the morphological variables of our 
studied lizards were not correlated with perch height. 
Thus, our ninth expectation is not fulfilled. On the oth- 
er hand, each population exhibited a unique pattern of 
relationship between thermal variables and perch height. 
These patterns are similar to those found in S. occidenta- 
lis, in which the body and environmental temperatures 
of some populations were independent of perch height 
(ADOLPH 1990), whereas in others, thermal characteris- 
tics were a primary factor modulating perch use (AsBURY 
& ADOLPH 2007). Our studied lizards of the GP used 
higher perches (up to two metres above ground) to behav- 
iourally lower their T, (especially around midday when 
the ambient temperature is higher on the ground), where- 
as lizards in the CP and TBP used higher perches (two 
metres above ground) to attain a higher T, (during the 


298 


morning when the environment is cooler on the ground; 
STEVENSON 1985). 

Results of interpopulational variation in perch height 
failed to confirm our tenth expectation and match a pattern 
found in Anolis aeneus, which was found to perch higher in 
a closed forest than it would in more open habitats, because 
of better access to insolation in the former that is not avail- 
able at lower heights in that habitat (Harris et al. 2004). 
As mentioned previously, each population selects a specific 
perch height according to their thermal environment (ex- 
cept in the TP), which directly affects their T,. Compara- 
tive and experimental studies in other Anolis species have 
found that lizard perch height also changes with the habi- 
tat in the presence of conspecifics and potential predators, 
triggering changes in behaviour and morphology, both of 
which improve its chances of survival (Losos et al. 1993, 
2004). Sceloporus grammicus shares its habitat with three 
lizard species and one snake species in the TP, whereas six 
lizards and two snakes co-occur in the TBP, just one lizard 
and two snakes in the CP, and two lizards and three snakes 
in the GP (Supplementary Table S7). Thus, it is likely that 
the intensity of competition and predation on S. grammi- 
cus would vary between the studied populations, and as a 
result, the lizards respond by perching at a height that al- 
lows them to avoid or diminish these biotic interactions 
as to enhance their chances of survival. Nonetheless, oth- 
er behaviours such as scanning for specific prey, searching 
for mates, and territoriality, could be determinants of the 
height at which S. grammicus perch in each population, as 
occurs in other lizard species (ANDREWS 1971). 

In conclusion, throughout its wide distribution range, 
S. grammicus encounters a great variety of combinations 
of environmental factors, suggesting that not all popula- 
tions experience the same degree of selective pressures. 
Therefore, we could not expect that different types of 
traits will diverge in the same trajectory. Our morpholog- 
ical comparisons demonstrated that the direction of sex- 
ual selection in males is similar in all populations, while 
the intensity of fecundity advantage is greater in females 
of the TP and GP than those in the CP and TBP. This latter 
selective pressure is of remarkable importance to popula- 
tion dynamics, which finally determines its own persis- 
tence over time and space. Our thermal ecology analy- 
ses demonstrate that intersexual variation is the excep- 
tion rather than the rule and, despite marked interpopula- 
tional differences in morphology, seasonality and habitat 
type play a more relevant role in the thermoregulation of 
this lizard species. In this regard, the plasticity of thermal 
traits in S. grammicus is complex and could be beneficial 
in the face of the current events of climate change and 
global warming. Perch height usage is not influenced by 
sex or by season, but exceptions do occur even if they are 
rare, as our results indicate. Overall, these results improve 
our knowledge of the ecology of a widely distributed liz- 
ard species, especially of the genus Sceloporus, of which 
many species continue to be little known to herpetologists 
despite their being quite diverse and abundant in North 
America. 
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